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ABSTRACT 

We present data from high-dispersion echelle spectra and simultaneous uvby photome- 
try for 7 Doradus. These data were obtained from several sites during 1994 November 
as part of the MUSICOS-94 campaign. The star has two closely-spaced periods of 
about 0.75 d and is the brightest member of a new class of variable early F-type stars. 
A previously suspected third period, very close to the other two, is confirmed. Previous 
observations indicated that sudden changes could be expected in the spectrum, but 
none were found during the campaign. The radial velocities rule out the possibility 
of a close companion. The phasing between the radial velocity and light curve of the 
strongest periodic component rules out the starspot model. The only viable mecha- 
nism for understanding the variability is nonradial pulsation. We used the method of 
moments to identify the modes of pulsation of the three periodic components. These 
appear to be sectorial retrograde modes with spherical harmonic degrees, {£,m), as 
follows: /i = (3,3), fa = (1,1) and fa = (1,1). The angle of inclination of the star is 
found to be i w 70°. 

Key words: Stars: pulsation - Stars: variables - Stars: individual: 7 Dor. 



1 INTRODUCTION 

7 Doradus (F0V) is the brightest member of what appears to 
be a new class of pulsating variable stars. In recent years, a 
number of late A- and early F-type dwarfs have been found 
to be variable with periods of the order of one day (see 
Krisciunas & Handler, 1995 for a list). This is too long for 
membership of the 8 Scuti class, where the periods are just 
a few hours. The period is, in many cases, compatible with 
that expected for rotation, which suggests that the variabil- 
ity may be due to the rotation of a spotted star (Mantegazza, 
Poretti & Zerbi 1994). However, most stars appear to be 
multiperiodic and require improbably large differential ro- 
tation. The starspot model for 7 Dor has been examined by 
Balona, Krisciunas & Cousins (1994b) who find that large 
overlapping spots would be required to give the observed 
light and colour amplitudes. 

The most likely explanation for the variability in these 



stars is nonradial pulsation (NRP), since it is difficult to 
understand multiple periodicities in terms of the starspot 
model or as orbital motion. Multiple periods are typical 
of NRP, but one needs to prove that the physical param- 
eters determined from an NRP model are realistic and that 
the other models are unable to explain the observations. In 
NRP, the light variability is due to the periodic variation 
in temperature and, to a smaller extent, the change in the 
shape of the star during pulsation. Pulsation with a period 
of about one day implies that the dominant restoring force 
is gravity (g-modes), rather than pressure (p-modes) as in 
the case of the S Set stars. Furthermore, because the stars 
show substantial light and radial velocity variations, the av- 
eraging effect over the visible hemisphere of the star must 
be quite small. This implies that the spherical harmonic de- 
gree, £, of the pulsation must be low, probably £ < 4. Aerts 
& Krisciunas (1996) have analysed 9 Aur, a star belonging to 
this group, using photometry and the cross-correlation pro- 
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file obtained with the CORAVEL instrument. They deduce 
that 9 Aur is pulsating in two modes, both having spherical 
harmonic indices (£, m) = (3, |1|). 

The mechanism which excites the pulsation is not 
known at present. The stars appear to be located at, or 
near, the red edge of the S Set instability strip. This in- 
dicates that the mechanism may be the same as for 5 Set 
stars (the ionization of H and He), but perhaps modified by 
convection. 

It would certainly be interesting to examine the line 
profile and light variations of these stars. These data may 
enable us to distinguish between NRP and rotational modu- 
lation of a spotted star. NRP produces line profile and light 
variations which, in principle, allow the determination of the 
spherical harmonic degree, I and the azimuthal number, m. 
This can be done by examining the variation of the moments 
of a line profile and calculating a discriminant. The value of 
the discriminant as a function of angle of inclination is used 
to determine the most probable mode {I, in). The method 
was developed by Balona (1986a, 1986b, 1987) and extended 
by Aerts (1993). 

The two periods of 7 Dor are stable and well determined 
from extensive photometric observations: Pi = 0.75701 d, 
P 2 = 0.73339 d (Balona et al. 1994a). More recently, evi- 
dence has been found for a third period at P4 = 0.67797 
d (Balona, Krisciunas & Cousins, 1994b). (We call this pe- 
riod P4 and not P3 because P3 was used in Balona et al. 
(1994a) to refer to a possible single period, later shown to 
be incorrect). Since the star is bright (V = 4.3) and well 
placed, it was chosen as one of the targets for the MUSICOS 
94 programme. The aim was to obtain a large number of 
high signal-to-noise line profiles which, in conjunction with 
Stromgren photometry, would decide which of the two mod- 
els, starspots or NRP, is the correct one and to determine 
the spot geometry or pulsation modes. Another factor in 
choosing to observe this star is the rapid change in the spec- 
trum which appears from time to time (Balona et al. 1994a). 
7 Dor is an IRAS source with a probable circumstellar dust 
cloud (Aumann 1985), which may be responsible for this 
effect. 

In this paper we present details of the observations and 
reduction procedure and analyse the data in terms of the 
two models. 



2 OBSERVATIONS 
2.1 Spectroscopy 

Spectroscopic data were obtained from four sites during 
1994 November: ESO (Chile), Mount Stromlo (Australia), 
Sutherland (South Africa) and Mount John (New Zealand). 
The projected rotational velocity of 7 Dor is high (v sin i — 
50 km s _1 ), which means that it is impossible to find a 
totally unblended spectral line. The most suitable line is 
probably that of Fe II at A5018.450 A or possibly Mg I 
at A5183.619 A. Instructions were given to ensure that the 
A5018.450 Fe II line was observed whenever possible. The 
time-coverage at the various sites is shown in Table 1. 

The spectra from ESO were obtained using the Coude 
Echelle Spectrometer (CES) fed by the 1.4-m CAT. The de- 
tector was CCD No. 9, a 1024 x 640 thinned back-illuminated 



Table 1. Details of the MUSICOS 94 campaign on 7 Dor. The 
epoch of the Julian day (JD) is 2449600.00; N is the number of 
spectra obtained. 



Site 


Observer 


JD start 


JD end 


N 


ESO 


Lagrange 


59.58 


59.88 


80 


ESO 


Lagrange 


60.52 


60.88 


113 


ESO 


Lagrange 


62.52 


62.87 


88 


ESO 


Foing 


63.58 


63.88 


90 


ESO 


Foing 


64.60 


64.88 


113 


ESO 


Foing 


65.61 


65.87 


79 


ESO 


Foing 


66.61 


66.90 


99 


MSSSO 


Lawson 


67.02 


67.09 


6 


MSSSO 


Lawson 


67.92 


68.05 


5 


MSSSO 


Lawson 


68.93 


69.04 


5 


MSSSO 


Lawson 


70.91 


71.13 


15 


MSSSO 


Lawson 


71.93 


71.98 


4 


MSSSO 


Lawson 


72.91 


73.13 


13 


MSSSO 


Lawson 


73.92 


73.92 


1 


MSSSO 


Lawson 


74.92 


75.19 


14 


Mt John 


Hcarnshaw 


72.92 


73.08 


6 


MtJohn 


Hcarnshaw 


74.90 


75.00 


5 


SAAO 


Bohm 


75.30 


75.60 


28 


SAAO 


Bohm 


76.31 


76.59 


25 


SAAO 


Bohm 


77.30 


77.60 


27 


SAAO 


Bohm 


78.28 


78.42 


11 


SAAO 


Bohm 


80.29 


80.31 


3 


SAAO 


Bohm 


81.28 


81.29 


2 


SAAO 


Bohm 


82.29 


82.30 


2 


SAAO 


Bohm 


83.29 


83.45 


5 


SAAO 


Bohm 


85.30 


85.30 


1 
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Figure 1. An example of the normalised spectrum of 7 Dor ob- 
tained at ESO. The feature at 5032 A is an artifact; the spectrum 
has not been smoothed. 

RCA-type CCD. The readout noise is about 25 e~ (2 pixel 
on-line horizontal binning) . The grating efficiency was about 
97 percent of maximum: the projected slit width was 48 /im. 
Only one order was observed with spectral coverage 4998 - 
5035 A and a sampling of 0.033 A pixel -1 . The exposure 
times were in the range 3-5 min, producing well-exposed 
spectra near the saturation limit. The resolving power was 
R — 40000. A Th-Ar lamp was used for wavelength cal- 
ibration. An example of the spectrum, normalised to the 
continuum (estimated by eye), is shown in Fig. 1. 

The Mt Stromlo and Siding Springs Observatories 
(MSSSO) spectra were obtained with the coude echelle on 
the 1.9-m telescope. The echelle grating has 79 lines mm" 1 
with the cross-disperser blazed at 5000 A. The detector was 
a 2048 x 2048 Tek CCD with gain of 2 e~/ADU and readout 
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Figure 2. An example of the normalised spectrum of 7 Dor ob- 
tained from several echelle orders at MSSSO. Spectra from SAAO 
look almost identical. A running mean box filter with a full width 
of seven pixels has been applied. 



noise of about 2600 e~ (no on-chip binning). The slit width 
was 200 fim (0.8 arcsec on the sky or 2 pixels on the detec- 
tor). The pixel size was about 25 /im and resolving power 
R w 60000. A Th-Ar lamp was used for wavelength cali- 
bration. The full spectral range of each order was covered. 
Each order covers about 87 A (0.043 A pixel" 1 ) and over- 
laps with neighbouring orders. An example of a spectrum 
obtained from 12 orders in the best-exposed range (4900 - 
5500 A) is shown in Fig. 2. 

The South African Astronomical Observatory (SAAO) 
spectra were obtained with the MUSICOS spectrograph at- 
tached to the 1.9-m telescope at Sutherland. Details of the 
instrumentation can be found in Baudrand & Bohm (1992). 
The detector was a 1024 x 1024 Tek CCD with gain of 2 
e~/ADU and readout noise of 46 e~ (no on-line binning). 
The pixel size was 24 jj,m; the spectral resolution R = 31000. 
The full spectral range of each order was covered, neighbour- 
ing orders overlap in wavelength. Each order covered about 
76 A (0.074 A pixel -1 ). The resulting spectra are very sim- 
ilar to those obtained as MSSSO (Fig. 2). A Th-Ar calibra- 
tion lamb was used. 

Spectra from Mt John University Observatory were ob- 
tained using the 1.0-m McLellan reflector and a Thomson 
TH 7882 CDA CCD as detector. This chip has 384 x 576 



pixels of dimension 23 fim. The slit width was 100 /jm, cor- 
responding to 2.0 arcsec on the sky and giving a wavelength 
resolution of about 0.2 A and a resolving power R — 25000. 
Details of the spectrograph are given in Hearnshaw (1977) 
and of the CCD detector system in Tobin (1992). These 
spectra have lower signal to noise ratio than the others and 
a running mean box filter with full width of 15 pixels was 
applied. Neighbouring orders do not overlap, so it is not pos- 
sible to construct a continuous spectrum as was done for the 
Mt Stromlo and SAAO data. The Mt John data were conse- 
quently used only to study the A5018.450 A line of Fe II and 
to monitor possible rapid changes in this and other spectral 
regions. 



2.2 Photometry 

Stromgren uvby photometry was obtained by L Balona 
during the last two weeks of 1994 November (simultane- 
ously with the MUSICOS run at Sutherland) using the 
single-channel photometer attached to the 0.5-m reflector 
at Sutherland. During this period, K Krisciunas obtained 
Johnson V photometry from CTIO using the 0.6-m Lowell 
reflector. The weather was very good from both sites: 12 
out of 14 nights were completely or partly photometric at 
Sutherland, while all 10 nights were photometric at CTIO. A 
total of 132 uvby and 245 V observations were obtained. The 
coverage was 10 - 18 observations per night at Sutherland 
and 20 - 30 observations per night at CTIO. In both cases 
a neutral density filter had to be used to avoid saturating 
the photomultiplier. The nearby F-type stars, HR 1291 and 
HR 1365 were used as local comparison stars. Observations 
of these two stars flanked each observation of 7 Dor. 

The spectroscopic and photometric data used in this 
paper can be obtained by examining the World Wide Web 
home page of the SAAO (www.saao.ac.za). 



3 SPECTROSCOPIC REDUCTIONS 

Reductions of the Mt John data were made using the MIDAS 
software package. All other data were reduced using SPEC1 
(single-order spectra) and SPEC2 (multi-order echelle spec- 
tra). These programs were written in FORTRAN by L Balona. 
The programs require, for each night, a list of files of three 
types: flatfield, arc and object frames. A few object frames 
are co-added and used to define the ridge line of the spec- 
trum as a function of (X, Y) position on the CCD. The 
flatfield frames can also be used for this purpose. A cer- 
tain number of pixels on either side of the ridge line are 
used to define the extent of the spectrum perpendicular to 
the dispersion axis (the X-axis). These were chosen so as 
to include practically all the light from the star. A one- 
dimensional spectrum is formed by adding the pixels along 
the Y direction inside the strip for a given X and subtract- 
ing the background for this value of X. The background is 
obtained from two strips on either side of the ridge line and 
sufficiently distant from it as to exclude any starlight. The 
two strips are added and a running mean applied to define 
the background as a function of X. 

The flatfield spectrum is obtained by co-adding all the 
flatfield frames (excluding those with saturated pixels) and 
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extracting a strip coincident with the ridge line of the stel- 
lar spectrum and of the same width. The pixels along the Y 
direction inside the strip are added and the background sub- 
tracted to form a one-dimensional spectrum. The resulting 
flatfield spectrum is normalised to unit mean and applied to 
the object spectra. 

The arc frames are co-added and a strip coincident with 
the ridge line of the stellar spectrum, and of the same width, 
is extracted. This is converted to a one-dimensional spec- 
trum and the background subtracted as described above. 
Finally, the flatfield correction is applied. The arc lines are 
fitted by Gaussians using non-linear least squares. This gives 
precise X positions, which, together with the laboratory 
wavelengths, are used to calibrate wavelength as a function 
of X. The arc lines are identified automatically by supplying 
the approximate wavelengths at two points on the spectrum 
as well as a list of wavelengths for the fitted lines. This infor- 
mation is contained in the configuration file. A polynomial 
of second-order was found sufficient to fit the calibration 
curve with a precision of 0.02 - 0.05 A rms. The output 
of SPEC1 is a list of ASCII files, one for each object frame, 
giving intensity versus wavelength. 

SPEC2 works in the same way as SPEC1, but on as many 
orders as desired. As input, the approximate Y positions 
of the centre of each order is required. Output is a binary 
file containing the intensity versus wavelength of each order. 
Subsequent programs read these data and calculate the re- 
quired intensity factors to match the orders (using the wave- 
length overlap between orders). Spikes are removed in the 
process to produce a continuous spectrum covering many or- 
ders. Finally, the continuum is estimated by eye and fitted to 
several points using a spline curve. This is used to normalise 
the spectrum. For consistency, the continuum points were 
chosen to be at the same wavelengths for all spectra. The 
Mt Stromlo spectra suffer from a misaligned flatfield, which 
produces a rather bumpy continuum for each order when 
the flatfield is applied. The normalisation technique just de- 
scribed was successful in correcting this problem. For both 
the Mt Stromlo and SAAO data we applied a box running 
mean filter with total width of seven pixels to remove excess 
noise (though the unsmoothed spectra could be used). 



4 PERIOD ANALYSIS 

The periodogram of the combined SAAO y and CTIO V 
photometry is shown in Fig. 3. The time span of 13 d is not 
sufficient to resolve the two closely-spaced frequencies. The 
single peak at / — 1.342 d _1 is a mean of fi = 1.32098 
and fi = 1.36354 d _1 . If /i is fitted to the data, a peri- 
odogram of the residuals has the strongest peak at 1.387 
d _1 . Conversely, if fi is fitted and removed from the data, 
the strongest peak in the periodogram is at 1.303 d _1 . Be- 
cause the total duration of the observations is too short to 
resolve f\ and fa , it is not possible to obtain accurate values 
of the frequencies. Nevertheless, the fact that prewhitening 
by the known value of /i produces a peak close to fi and 
vice- versa, shows that the two frequencies are still present. 

Since /i and fi are still present, it is reasonable to as- 
sume that the frequencies have remained unchanged over the 
last few years. To test this, we show in Table 2 the ampli- 
tudes and phases for individual seasons using /i = 1.32098 
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Figure 3. Periodogram of the combined SAAO y and CTIO V 
data of 7 Dor (top panel) and the periodogram after removal of 
the two principal components, fi and fi (bottom panel). Ampli- 
tude is in mmag, frequency in cycles d _1 . 

and fi — 1.36354 d -1 (the third frequency, fi, is discussed 
below). It can be seen that the phases over the five years do 
not differ significantly from the mean, in spite of the small 
formal errors (these do not take into account the distribu- 
tion of the observations). The maximum deviation in phase 
is 0.09 periods for /i and 0.06 periods for /2, while the stan- 
dard deviation is about 0.03 periods in both cases. The low 
amplitude during the 1993 season is particularly striking. 

Fig. 3 also shows the periodogram of the 1994 Novem- 
ber y and V data after removing the two periodic compo- 
nents. It is clear that a third frequency at / = 1.4888 d _1 is 
present. This was strongly suspected by Balona, Krisciunas 
& Cousins (1994b) who obtained / 4 = 1.475 d"\ A A = 7.4 
mmag. This third frequency, f±, is therefore confirmed. Ta- 
ble 2 shows that the amplitude of this component was below 
the detectable threshold (about 3 mmags) prior to the 1993. 
A periodogram of the combined 1993 and 1994 seasons gives 
a best estimate of / 4 = 1.47447(±5) d _1 . Table 3 lists the 
amplitudes and phases of the three periodic components of 
the V and Stromgren colour indices for the combined 1993 
and 1994 seasons. This gives the best estimate of the phases 
for 1994. However, the amplitudes of /i and fi for 1994 
should probably be taken from Table 2, since they were par- 
ticularly low in 1993. 



5 RADIAL VELOCITY 

The radial velocity is determined from the displacement of 
a spectral line from its laboratory wavelength. If the line 
is symmetric, the displacement is uniquely defined, since 
the centroid, median and mode coincide. If it is not sym- 
metric, problems arise as to what is used in defining the 
displacement. This is always a problem in pulsating stars, 
since the lines are almost always skew. The cross-correlation 
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Table 2. Amplitudes (Aj) and phases (<pj) of a three-component Fourier fit to the V data between 1989 and 1994. The frequencies are 
fi = 1.32098, Si = 1-36354, / 4 = 1.47447 d" 1 and the epoch is JD 2449000. Semi-amplitudes arc in millimags; phases are in fractions 
of a period. The JD is with respect to 2440000; N is the number of observations. 





JD 


N 


Ai 


01 


A 2 


2 


Ai 


<f>4 


7838 - 


7970 


175 


16.8 ± 1.0 


0.03 ±0.01 


12.9 ± 1.0 


0.26 ± 0.02 


3.0 ±1.0 


0.56 ± 0.05 


8185 - 


8327 


107 


12.8 ±1.2 


0.04 ±0.02 


15.6 ± 1.2 


0.27 ±0.01 


2.6 ±1.2 


0.39 ± 0.08 


8993 - 


9005 


278 


16.3 ±1.2 


0.09 ±0.01 


19.5 ± 1.5 


0.31 ±0.01 


2.6 ±0.3 


0.26 ± 0.03 


9326 - 


9406 


331 


6.8 ±0.7 


0.98 ±0.02 


9.5 ±0.7 


0.19 ±0.01 


5.9 ±0.6 


0.91 ± 0.02 


9672 - 


9685 


377 


14.9 ±0.7 


0.92 ±0.01 


13.4 ±0.7 


0.21 ±0.08 


7.2 ±0.5 


0.14 ±0.01 



Table 3. Amplitudes (Aj) and phases (<j>j) of a three-component Fourier fit to the radial velocity, V r and V, b — y, ci and u — b (1993 
and 1994 seasons). The frequencies are /i = 1.32098, / 2 = 1.36354, / 4 = 1.47447 d~ x and the epoch is JD 2449000. Semi-amplitudes 
are in km s _1 for V r and in millimags for the photometry; phases are in fractions of a period. 







01 


A 2 


2 


A 4 


04 


V r 


0.3 ±0.1 


0.02 ± 0.05 


1.3 ±0.1 


0.64 ±0.01 


0.6 ±0.1 


0.47 ± 0.02 


V 


11.6 ±0.5 


0.97 ±0.01 


13.5 ±0.5 


0.21 ±0.01 


6.7 ±0.5 


0.14 ±0.01 


b~y 


3.9 ±0.5 


0.02 ± 0.02 


2.5 ±0.5 


0.32 ±0.03 


2.1 ± 0.4 


0.21 ±0.03 


Cl 


11.7 ±0.8 


0.53 ± 0.01 


7.9 ±0.8 


0.80 ±0.02 


2.5 ± 0.7 


0.63 ±0.05 


u — b 


5.8 ±0.7 


0.55 ± 0.02 


4.1 ±0.7 


0.85 ±0.03 


1.7 ±0.6 


0.60 ±0.06 



technique is often used to determine the radial velocity, but 
this does not get round the problem which is now trans- 
ferred to defining the position of the (possibly asymmetric) 
correlation function. In this paper we will adopt the centroid 
(the first moment of the line profile) as the measure of dis- 
placement because it is easy to calculate and because it is 
used in theoretical work on pulsating stars. 

In the case of 7 Dor, the most accurate measurement 
of radial velocity is obtained by making use of the echelle 
spectra from MSSSO and SAAO. We used the first spec- 
trum obtained at SAAO as the master and determined the 
radial velocity by cross-correlation with this spectrum. The 
periodogram of the radial velocities (Fig. 4) shows a peak at 
/ = 1.365 ± 0.002 d"\ very close to f 2 . Although the data 
sampling is not sufficient to resolve the three frequencies 
in the periodogram, numerical tests show that the correct 
amplitudes are recovered if the frequencies are given. Since 
we already know the frequencies with high accuracy, Fourier 
decomposition of the radial velocities allows the amplitudes 
and phases to be determined: these are shown in Table 3. 
The light, colour and radial velocities for each of the three 
periodic components are shown in Fig. 5. The standard de- 
viation of one observation for the radial velocity is found to 
be 0.68 km s" 1 . 

An important result of these observations are the ratios 
of radial velocity to light amplitude for the three periodic 
components. Note that the light amplitudes for the /1 and 
f'2 components are comparable, but that the radial velocity 
amplitude of /1 is very much smaller than that of / 2 . It 
is not easy to understand how two starspots could produce 
such different radial velocity to light amplitude ratios. On 
the other hand, NRP can account for this difference if the 
two components belong to different modes. The ratios of 
radial velocity to light amplitude for / 2 and / 4 are nearly 
the same, indicating that (£, m) for these two components 
could be the same (as indeed they are - see below). Note 
also that in / 2 , radial velocity maximum occurs near light 
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Figure 4. Periodogram of the radial velocities obtained by cross- 
correlation of the SAAO and MSSSO multi-order echelle data (top 
panel) and after removing the frequency / = 1.365 d _1 (bottom 
panel). The amplitude is in km s , the frequency in cycles d . 



maximum. This again is not consistent with the starspot 
model, where you would expect radial velocity maximum to 
occur when the spot is receding, i.e. a quarter of a period 
after minimum light for a dark spot (or a quarter of a period 
after maximum light for a bright spot). 
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Figure 5. The radial velocity, light and colour curves for the periodic component /i = 1.32098 d _1 (left), /2 = 1.36354 d _1 (middle) 
and /4 = 1.47447 d _1 (right). The epoch of phase zero is JD 2449000. Tick marks arc spaced at intervals of 2 km s~ 1 for the radial 
velocity curve and 0.02 mag for the light and colour curves. 



6 THE MOMENTS 

To determine the pulsation modes for each of the three pe- 
riods in 7 Dor, we need to measure the first few moments of 
a particular spectral line as a function of time. As already 
mentioned, the A5018.450 A line of Fe II was chosen be- 
cause it is strong and relatively unblended. The procedure 
for obtaining the moments has been automated in a program 
called DoSpec (written by L Balona). 

In general, DoSpec looks for spectral absorption lines 
with central depth greater than a certain threshold, fits them 
by Gaussians, and removes them from the spectrum. The 
threshold is lowered and the procedure repeated until a cer- 
tain level is reached. DoSpec then adds each fitted line back 
to the spectrum, calculates the moments, and removes the 
line again. In this way moments are calculated which are 
relatively free of the effects of neighbouring lines. The out- 
put is a file listing all lines in the spectrum and their first 
seven moments (moments 2-6 are calculated relative to the 
centroid). The moment of order zero is just the equivalent 
width (EW), the first moment is the centroid (from which 
the radial velocity can be calculated), the second moment is 
the rms width of the line, etc. 

DoSpec was applied to the neighbourhood of the A5018 
line using all available spectra (ESO, Mt Stromlo, SAAO 
and Mt John). The periodogram of the radial velocity of 
this line from the combined data is rather puzzling as the 
strongest peak is at / = 0.734 d _1 . This is different from 
any frequency so far identified, though it is rather close to 
/4/2. However, if only the SAAO and Mt Stromlo data are 
used, the periodogram is quite different, having a peak at 
/ = 1.359 d _1 . This is clearly /2, identified above from 
the radial velocity obtained by cross-correlation. The pe- 
riodogram from the ESO and Mt John data alone shows a 
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Figure 6. Periodogram of the first moment of the A5018.450 A 
line of Fe II. Top panel - all available data; middle panel - MSSSO 
and SAAO cchcllc data; bottom panel - ESO and Mt John data. 
The amplitude is in km s _1 , the frequency in cycles d _ . 



peak at / = 0.728 d 1 . These periodograms are shown in 
Fig. 6. 

We do not understand why the ESO data do not show 
any power at or near / = 1.36 d~ .A phase diagram of 
the first moment (mi) of the ESO data phased with /i, 
and /4 (Fig. 7) shows reasonable agreement with the SAAO 
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Figure 7. Phase variation of the first moment of the Fe II line 
A5018.450 A from the ESO data for the three periodic compo- 
nents. The epoch of phase zero is JD 2449000. Tick marks are 
spaced at intervals of km s _1 . 

& MSSSO curves (see Fig. 8). But for /2 the curve seems 
to be of the double-wave type with a minimum at about 
phase 0.9 and another minimum at phase 0.4, where the 
SAAO & MSSSO data show a maximum. The minimum 
at phase 0.4 arises almost entirely from data for one night 
(JD 2449664). There is no reason to suspect an instrumental 
cause for this behaviour, but unfortunately we do not have 
photometry taken at about the time when the double-wave 
form developed for fa. The star is known to display irregular 
excursions in light from time to time (Balona et al. 1994a). 
It is possible that such an event may be responsible for the 
minimum where a maximum is expected, but this cannot be 
confirmed. One should also take into account the fact the 
ESO data are rather poorly sampled, several nights having 
similar mean phase for /i or /2 . In view of this problem, we 
decided it would be better to use the SAAO & MSSSO data 
only in determining the modes. 



7 MODE IDENTIFICATION 

Since the phase relationship between the light and radial 
velocity variations rules out the starspot hypothesis, and 
since we cannot think of any model invoking orbital motion 
which may account for the variations, it is safe to assume 
that NRP is the most probable explanation. However, we 
need to prove that the parameters required for NRP are 
physically realistic. To do this, we need to identify the modes 
of oscillation and determine these parameters. 

Any arbitrary distribution can be represented by speci- 
fying its moments. The greater the number of moments, the 
more accurate is the representation. During pulsation, the 
line profile will vary and so will its moments. For example, 
in 7 Dor the moments will vary with the three frequencies, 



/i, /2 and fi. The amplitude and phase of a particular mo- 
ment which belongs to one of the three oscillations can be 
calculated by fitting a multiperiodic Fourier curve. For a 
particular angle of inclination, i, and a given mode, (£, m), 
the amplitudes and phases of the first few moments can be 
fitted to the predicted ones, and the pulsation parameters 
obtained by least squares. In general, the fit will be poor 
unless the chosen values of i, £ and m are close to the true 
values. In this way, the most likely values of i, £ and m can 
be estimated. 

This method of mode estimation was developed by 
Balona (1986a, 1986b, 1987). It uses the algorithms de- 
scribed in Balona (1987) and assumes that the pulsation 
velocity amplitude is much smaller than the projected ro- 
tational velocity. This is certainly true in 7 Dor where the 
radial velocity semi-amplitude is less the 2 km s _1 , whereas 
vsmi = 50 km s^ 1 . The algorithm uses a first-order ex- 
pansion in the ratio fi/w (SI being the frequency of rotation 
and ui the frequency of pulsation). The eigenfunctions of a 
rotating, pulsating star are well described by pure spherical 
harmonics only when this ratio is small. This is not the case 
in 7 Dor. In Balona et al. (1994a), the stellar radius is found 
to be R = 1.32iio. Using the value v sin i — 50 km s _1 , gives 
a rotation period P ro t < 1-34 d. Therefore SI/lo > 0.55. For 
such a large ratio, departures of the eigenfunctions from pure 
spherical harmonics are bound to be significant and we do 
not expect a very good fit on the basis of a first-order ex- 
pansion. Spheroidal modes with different values of (£, m) are 
mixed in, and so are toroidal modes. These additional com- 
ponents will have lower amplitudes, so we can expect the 
method to identify the dominant spherical harmonic mode 
(£,m). 

Aerts & Krisciunas (1996) apply a similar method in 
their analysis of 9 Aur. However, the rotational velocity of 
this star is rather low and the velocity variations larger than 
in 7 Dor. The approximation discussed above is probably 
not a good one for this star, and a different algorithm was 
used. Owing to the low rotational velocity, the algorithm is 
incapable of distinguishing between prograde (m < 0) and 
retrograde (m > 0) modes. 

We use the Fe II line A5018.450 A to probe the pulsa- 
tion mode. It may be possible to use more than one line to 
improve the signal to noise ratio. In fact, Aerts & Krisciunas 
(1996) use the profile obtained from the cross-correlation of 
the spectrum with a mask in the CORAVEL instrument for 
this purpose. We have not tested the reliability of this pro- 
cedure. First, we calculate the centroid (first moment) of the 
line profile, which gives the radial velocity. Then we calcu- 
late higher-order moments with respect to the centroid as 
origin. The moments which are actually required are those 
with respect to the centre of mass of the star. This means 
that we need to determine the mean radial velocity, and 
then re-calculate all the higher order moments with respect 
to this velocity (rather than the centroid) . Transformation of 
the origin of the moments from the centroid to the centre-of- 
mass velocity is easily accomplished using a polynomial ex- 
pansion. Using only the MSSSO and SAAO data, the mean 
radial velocity is < mi >= 22.4 ± 0.2 km s _1 ; the standard 
deviation of one observation being 2.3 km s _1 . The standard 
deviation is larger than in the cross-correlation technique 
because only one line is used. A few published radial veloc- 
ity measurements are available (Abt & Biggs 1972). These 
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range from 20 to 27 km s~\ in good agreement with the 
value obtained here from the Fe II line. 

In Fig. 8, we show the variation of the equivalent width 
(EW) and the first four moments of the Fe II line A5018.450 
A phased with /i, f-2 and fi for the combined MSSSO and 
SAAO data. The amplitudes and phases are shown in Ta- 
ble 4. We have fitted the data with all three periodicities 
because these periodicities are present in the light curve. 
It might be argued that the only period which is definitely 
present in the radial velocities is /2. It is certainly true that 
/2 has the highest amplitude, but this could be merely a 
result of the particular nonradial mode favouring the first 
moment as opposed to other moments. The amplitude of the 
radial velocities (and other moments) places a constraint on 
the mode identification even if it is close to zero. It is there- 
fore necessary to determine the amplitudes and phases of the 
moments of all three periodicities if one is to understand how 
the light variations are produced. 

The mean equivalent width is < EW >= 0.369 ± 0.001 
A (the standard deviation for one observation is 0.015 A). It 
is surprising to find that there is significant variation of the 
EW when phased with /i, but not with fi or Ja- For the /i 
component, the EW reaches a maximum at about the time 
of maximum ci colour, i.e. slightly before maximum light. 
On the other hand, there is no significant variation with /2 
even though the light and colour amplitude is almost the 
same as /i. 

Only a few moments show significant phase variation, 
which is not surprising in view of the low radial velocity am- 
plitude of the star. This means that mode identification will 
be indeterminate unless we use the most significant varia- 
tions and employ as many constraints as possible. In the 
case of 7 Dor we need to consider how we can constrain the 
pulsational parameters. NRP generates a variable tempera- 
ture distribution on the photosphere which distorts the line 
profile owing to rotation of the star. This effect can be quan- 
tified as a fictitious velocity field, Vf. It is fictitious because it 
is not generated by the variation of displacement during pul- 
sation. Instead, the spatial variation of temperature across 
the photosphere caused by nonradial pulsation, in combina- 
tion with rotation of the star, adds or removes flux from the 
line profile. The net result is a change in the radial veloc- 
ity as measured by the centroid or line minimum. It plays 
an important role only if rotation is reasonably rapid. Very 
crudely, the amplitude of the velocity generated by this ef- 
fect is approximately given by the product of the relative 
luminosity variation and the projected rotational velocity. 
In 7 Dor, the semi-amplitude of the light variation is about 
0.01 mag and usini = 50 km s _1 , so we expect Vf ~ 0.5 
km s _1 . This is comparable to the observed radial velocity 
semi-amplitude and cannot be neglected. 

The other velocity fields are v r - the pulsational veloc- 
ity in the vertical direction and Vh - the pulsational velocity 
in the horizontal direction. For each of these three velocity 
fields, the amplitude and phase needs to be determined, giv- 
ing six free parameters for a given i, £ and m. We know that 
v r must be small compared to Vh because the three modes 
must all be g-modes. The ratio 

vh „ 4n 2 fR A 

v r ~ GM ' 



is approximately 36 for 7 Dor. Hence we may safely ignore 
v r and limit the solutions to Vf and Vh only, resulting in 
four free parameters. A further constraint can be imposed: 
Vf < v sin i, since larger values imply a relative flux variation 
greater than unity. 

To first order, the horizontal velocity component, v h , 
will not contribute to the light variation since there is no 
geometric distortion. All the light variation arises from the 
temperature variation during pulsation (vj). On the other 
hand, Vh will certainly contribute to the distortion of the 
line profile. The amplitude and phase of the light curve is 
determined only by Vf , but the amplitude and phase of the 
radial velocity curve depends on vh as well as Vf . This allows 
for a different phase relationship between the two curves 
other than that predicted by the spot model. 

The light curve and the radial velocity curve can there- 
fore be used to determine Vf and vn for any i, £ and m, 
but at least one more moment is required to produce a dis- 
criminant. In principle, a solution can be obtained using the 
amplitudes and phases of the light curve and the first two 
moments only. However, there is information in the third 
and higher moments which should be used to constrain the 
solution. It is clear, however, that the higher order moments 
should be given less weight because they have larger uncer- 
tainties. The weights to be assigned to a particular moment 
can be determined from the standard error of its amplitude 
and phase (easily calculated from the Fourier fit) . The equa- 
tions of condition in the algorithm have dimensions of ve- 
locity, since each equation is divided by (usini)™ -1 for the 
n-th moment. Hence the appropriate weight is obtained by 
dividing the standard error of the moment by this quan- 
tity. This results in the following standard deviations being 
assigned to each moment: light curve - 1.00; first moment 
- 2.09; second moment - 2.73; third moment - 6.82; forth 
moment - 11.84 km s _1 . 

The zero points of the second and fourth moment of 
the line profile are sufficient to determine the projected ro- 
tational velocity and the rms line width (Wi) of the intrinsic 
profile. Wi includes a contribution to the intrinsic line width 
owing to the finite resolution of the instruments (the differ- 
ence in resolving power of the SAAO and MSSSO spectro- 
graphs is equivalent to about 5 km s— 1). In principle, it is 
possible to determine Wi for the two instruments, but it is 
unlikely that this difference would be significant. The zero 
points of the second and fourth moments are rri2 = 1897± 14 
and 7714 = (10.54 ± 0.16) x 10 6 in units of km s _1 to the re- 
spective powers. Using a limb darkening value of u = 0.5, 
appropriate for an early F-type star in the visible region, we 
find vsini = 48 km s _1 and Wi = 37 km s _1 . The value of 
usini agrees well with that of Slettebak et al. (1975). These 
values are used to deconvolve the moments to the values 
they would have had in the limit of zero intrinsic width as 
described in Balona (1986a, 1986b, 1987). 

Using these values, we calculated the discriminant, a, 
for a choice of mode (I, m) as a function of angle of incli- 
nation, i, using the weighted light curve and first four mo- 
ments. In this calculation, we adopted a radius R = 1.32Rq, 
a limb darkening coefficient u = 0.50, and the Ledoux ro- 
tational splitting coefficient, C = 0.1 (a typical value; the 
exact value is unimportant). The computation included all 
modes up to £ = 4. The discriminant is plotted against angle 
of inclination in Fig. 9 for the three modes. 
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Figure 8. Phase variation of the equivalent width (EW) and the first four moments (mo ~ m i) 

of the Fe II line A5018.450 A phased 

with /i (left), f 2 (middle) and /4 (right). The epoch of phase zero is JD 2449000. Tick marks are spaced at intervals of 0.02 A for EW 
and 2 km s _1 for mi. The other moments are plotted on an arbitrary scale. 



Table 4. Amplitudes (Aj) and phases (<pj) of a three-component Fourier fit to the EW and first four moments of the Fc II line 
A5018.450 A for the combined MSSSO and SAAO data. The frequencies arc /i = 1.32098, f 2 = 1.36354, / 4 = 1.47447 d" 1 and the 
epoch is JD 2449000. Scmi-amplitudcs arc in mA for the EW, in km s~ x for mi and in powers of km s — 1 for the other moments 
(m2 X 10 2 ,m3 X 10 4 ,m4 X 10 5 ). Phases are in periods. 







<Pi 


A 2 


<i>2 


A 4 


4>4 


EW 


11.4 ±2.2 


0.60 ±0.03 


5.9 ±2.6 


0.24 ± 0.06 


5.0 ± 1.8 


0.04 ± 0.06 


mi 


1.0 ±0.3 


0.30 ±0.06 


2.3 ±0.4 


0.61 ± 0.03 


1.4 ±0.3 


0.53 ±0.03 


m-2 


1.1 ±0.3 


0.59 ±0.04 


0.4 ±0.3 


0.38 ± 0.12 


0.6 ±0.2 


0.12 ±0.05 


mz 


0.7 ±0.3 


0.27 ±0.07 


1.9 ±0.3 


0.62 ± 0.03 


1.2 ± 0.3 


0.54 ±0.03 


ni4 


9.0 ±3.0 


0.58 ±0.05 


4.8 ±3.0 


0.39 ± 0.11 


6.1 ± 2.4 


0.19 ±0.06 



Since the angle of inclination must be the same, f 2 con- 
strains i > 50°. The (3, 1) mode leads to a light amplitude 
which is too low except over a very small range near i = 77° . 
Therefore, we can identify f 2 as either (2, -1) or (1, 1). Both 
give realistic pulsation parameters and good agreement with 
the light amplitude for an inclination 65 < i < 80°, but the 
(2, -1) mode does not match the phase of the light curve very 
well. With this limit on the angle of inclination, we find that 
fi can only be (3, 3). This also gives realistic pulsation pa- 
rameters and good agreement with the light curve. Finally, 
we find (1, 1) as the best estimate for /4, since (2, 0) and 
(2, -1) produce light amplitudes which are twice as high as 
observed. 

Another way of imposing the constraint that the incli- 
nation angle for all three modes must be the same, is to 
calculate the sum of the weighted discriminant at a partic- 
ular value of i for the three modes for all values of (£,m). 
When the minimum of this sum is plotted against angle of 
inclination, an indication of the uniqueness of the solution 



can be ascertained. This is shown in Fig. 10 where we have 
chosen weights of unity for /i and f 2 , and a weight of 0.25 
for the less significant f^. There is no solution for i < 7° 
because this would imply an equatorial velocity larger than 
the breakup velocity. The variance for 45 < i < 63° is too 
large to be plotted. The figure shows that the best solution 
is found for i > 63°. As mentioned above, the identifications 
for this range are (3, 3) for /i, (2, -1) or (1, 1) for f 2 and 
(1, 1) for / 4 . 

On the basis of these data, we conclude that the angle 
of inclination of 7 Dor is probably i fa 70° , and that the best 
estimate of the values of (i, m) for the three oscillations are: 
(3, 3) for /1, (1, 1) for f 2 and for f^. Table 5 gives the 
pulsational parameters for this angle of inclination. 
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Table 5. Pulsation parameters for the most probable mode identifications, (£,m), for the three oscillations in 7 Dor. The angle of 
inclination is assumed to be i = 70°. 
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Figure 9. Discriminant (arbitrary units) as a function of angle 
of inclination (in degrees) for various modes (£,m). Top panel: 
/l = 1.32098 d _1 ; middle panel: f 2 = 1.36354 d _1 ; bottom panel: 
f 4 = 1.47447 d" 1 . 
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Figure 10. The minimum sum of the discriminant variance over 
all three modes is shown as a function of angle of inclination (in 
degrees). The best solution is found for 63 < i < 80°. 

8 DISCUSSION 

One of the most puzzling aspects of 7 Dor is the sudden ap- 
pearance, from time to time, of isolated absorption or emis- 
sion features (Balona et al. 1994a). Careful inspection of 
all the spectra obtained during the MUSICOS-94 campaign 
has failed to reveal such oddities. The star appears to have 
a normal early-F spectrum, though broadened by rotation. 



There is no doubt, however, that the light variations of the 
star cannot be explained entirely by three periodic compo- 
nents. There is a certain level of non-periodicity which is 
not yet understood (Balona et al. 1994a, Balona, Krisciu- 
nas & Cousins 1994b). The spectroscopic observations dur- 
ing MUSICOS-94 show evidence for this effect on one night 
(JD 2449664), where the radial velocities are substantially 
smaller than expected. 

Our new photometric observations confirm the exis- 
tence of a third oscillation, f 4 — 1.47447 d~ x , previously 
suspected by Balona, Krisciunas & Cousins (1994b). The 
existence of three frequencies implies substantial differen- 
tial rotation on the basis of the starspot model. This model 
also fails to account for the great differences which exist be- 
tween the radial velocity to light amplitude ratios between 
the two main oscillations. But the most serious objection 
to the starspot model is the fact that for f 2 radial velocity 
and light maximum coincide in phase. One would expect a 
90° phase difference for a starspot. The small radial velocity 
amplitude also excludes the possibility of a binary compan- 
ion. We conclude that nonradial pulsation is the only viable 
explanation for the observations. 

From the light curve and first four moments of the line 
profile of Fe II A5018.450 A, we show that 7 Dor has a 
high angle of inclination, probably i fa 70° and that the 
three modes may be identified as /1 = (3, 3), f 2 = f± = 
(1, 1) with fair confidence. The pulsation parameters indi- 
cate quite a large value Vf = 15.2 km s _1 for /1. Now the 
relative flux amplitude 

AF _ v f 
F v sin i 

which gives AF/F = 0.32 for /1, implying that this mode 
produces a large temperature variation AT /T fa 0.08. This 
may account for the significant phase dependence of the 
equivalent width for this mode, whereas the other modes 
show little or no phase dependence. The fact that the am- 
plitudes of the colour indices ci and (u — b) are larger in this 
mode is also consistent with the large Vf. 

The fact that f 2 dominates the radial velocity varia- 
tion, but that fi does not, though they have similar light 
amplitudes, is due to the fact that /1 has a relatively high 
spherical harmonic degree I = 3, whereas f 2 has the lower 
value 1=1. The averaging effect in the latter case is much 
less, leading to a larger radial velocity amplitude. On the 
other hand, the much larger value of Vf for /1 compensates 
in the light curve. 

It is interesting that all three modes seem to be sec- 
torial retrograde modes. The frequency of pulsation in the 
co- rotating frame, vo, is related to the frequency of pulsa- 
tion in the observer's frame, v, by: i>o = v + rnQ, where O 
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is the frequency of rotation. If i — 70° and R = 1.32i?Q, we 
find fi = 1.04 d _1 (P rot = 0.96 d). Thus, in the co-rotating 
frame, the three oscillations have the following frequencies: 
/io « 4.44, / 20 « 2.40 and / 40 ~ 2.51 d" 1 . /i and f 2 are 
substantially different and must have quite different radial 
orders in spite of the fact that they are nearly equal in the in- 
ertial frame. On the other hand, /2 and /<t, which are of the 
same degree, £, probably have neighbouring radial orders. 
This information could be used to extract seismological in- 
formation, but the effect of rotation needs to be taken into 
account more accurately before this can done. 

The success of mode identification in this star is encour- 
aging. The application of the method to other stars might 
well result in very interesting results, as in the case of 9 Aur 
(Aerts & Krisciunas 1996), but other types of pulsating star 
may be equally interesting. For example, an intensive cam- 
paign on a periodic Be star (and there are many of these 
which are bright) should result in a definitive test of the 
NRP/starspot models (Balona 1996). All that would be re- 
quired is an accurate determination of the radial velocity to 
light amplitude ratio. 
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